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, Abstract. The Maxwell-Boltzmannian approach to nuclear reaction rate theory is extended 

to cover Tsallis statistics (Tsallis, 1988) and more general cases of distribution functions. An 
analytical study of respective thermonuclear functions is being conducted with the help of sta- 
tistical techniques. The pathway model, recently introduced by Mathai (2005), is utilized for 
thermonuclear functions and closed-form representations are obtained in terms of H- functions 
and G-functions. Maxwell-Boltzmannian thermonuclear functions become particular cases of 
the extended thermonuclear functions. A brief review on the development of the theory of 
analytic representations of nuclear reaction rates is given. 

> . 

2^ ■ 1 Introduction 

In the evolution of the Universe, chemical elements are created in cosmological and stel- 
qq , lar nucleosynthesis (Clayton, 1983; Fowler, 1984). Solar nuclear energy generation and solar 

neutrino emission are governed by chains of nuclear reactions in the gravitationally stabi- 
lized solar fusion reactor (Davis, 2003). One of the first utilization of Gamow's quantum 
mechanical theory of potential barrier penetration to other than the analysis of radioactive 
decay was to the question on how do stars generate energy (Critchfield, 1972; for a brief 
essay on the history of this discovery see Mathai and Haubold, 1988). Continued attempts 
are aiming at generating energy through controlled thermonuclear fusion in the laboratory. 
In a nuclear plasma, the rate of reactions and thus energy release can be determined by an 
average of the Gamow penetration factor over the distribution of velocities of the particles 
of the plasma (Haubold and Mathai, 1984). Understanding the mathematical and statistical 
methods for the evaluation of thermonuclear reaction rates is one of the goals of research in the 
field of stellar and cosmological nucleosynthesis. Practically all applications of fusion plasmas 
are controlled in some way or another by the theory of thermonuclear reaction rates under 
specific circumstances. After several decades of effort, a systematic and complete theory of 
thermonuclear reaction rates has been developed (Haubold and John, 1978; Anderson et al., 
1994; Haubold and Mathai, 1998; Mathai and Haubold, 2002). Reactions between individual 
particles produce a randomization of the energy and velocity distributions of particles. The 
depletion of particles by reactions is balanced by the diffusion of the particles in the macro- 
scopically inhomogeneous medium. As a result of this balance, the fusion plasma may reach 
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a quasi-stationary state close to equilibrium, in which steady fluxes of matter, energy, and 
momentum are present. This also lead to the assumption that the distribution of particles 
can be assumed to be Maxwell-Boltzmannian in almost all cases of interest to stellar physics 
and cosmology. 

The derivation of closed-form representations of nuclear reaction rates and useful approx- 
imations of them are based on statistical distribution theory and the theory of generalized 
special functions, mainly in the categories of Meijer's G-function and Fox's H-function (Mathai 
and Saxena, 1973; Mathai, 1993; Aslam Chaudhry and Zubair, 2002). For an overview on 
the application and historical background of integrals and distribution functions for reaction 
rates and their representation in terms of special functions, see Hegyi (1999) and Moll (2002). 
Special cases which can be derived from the general theory through expansion of respective 
physical parameters (like the cross section factor) are nonresonant (Hussein and Pato, 1997; 
Ueda et al., 2000) and resonant (Ueda et al., 2004; Newton et al., 2007) reaction rates, reaction 
rates for cosmological nucleosynthesis (Bergstroem et al., 1999) and fitting of experimental 
data to analytic representations (Brown and Jarmie, 1990). Specific mathematical methods 
for deriving approximate analytic representations of nuclear reaction rates are expansions 
of hyper geometric functions (Saigo and Saxena, 1998), transformation of extended gamma 
functions (Aslam Chaudhry, 1999), and asymptotic expansion of the Laplace transform of 
functions (Ferreira and Lopez, 2004). 

Only recently, related to the production of neutrinos in the gravitationally stabilized solar 
fusion reactor, the question of possible deviations of the velocity distribution of plasma par- 
ticles from the Maxwell-Boltzmannian case due to memory effects and long-range forces has 
been raised (Coraddu et al., 1999; Lavagno and Quarati, 2002; Coraddu et al., 2003; Lissia 
and Quarati, 2005; Lavagno and Quarati, 2006). This was initiated by Tsallis' non-additive 
generalization of Boltzmann-Gibbs statistical mechanics which generates q-exponential func- 
tion as the fundamental distribution function instead of the Maxwell-Boltzmann distribution 
function. Tsallis statistics covers Boltzmann-Gibbs statistics for the case q = 1 (Tsallis, 1988; 
Gell-Mann and Tsallis, 2004; Tsallis, 2004). A first attempt was made to extend the theory 
of nuclear reaction rates from Maxwell-Boltzmann theory to Tsallis theory (Saxena et al., 
2004; Mathai, 2005; Mathai and Haubold, 2007). This paper develops the complete theory 
for closed-form representations of nuclear reaction rates for Tsallis statistics. In this context 
an interesting discovery was made by Mathai (2005), namely that even more general distri- 
bution functions can be incorporated in the theory of nuclear reaction rates by appealing to 
entropic and distributional pathways. Starting from a generalized entropy of order a, through 
the maximum entropy principle, distribution functions are generated which include Maxwell- 
Boltzmann and Tsallis distributions as special cases. 

In Subsections 1.1, 1.2, 1.3, and 1.4, we introduce the definition of the thermonuclear 
reaction rate and respective thermonuclear functions for the standard, cut-off, depleted, and 
screened case, respectively. Each subsection provides the integral of the thermonuclear func- 
tion for the cases of Maxwell-Boltzmann distribution and a distribution, the latter covers the 
q-exponential of Tsallis. Section 2 provides prerequisites for the use of G- and H- functions in 
Mellin-Barnes integral representation and also discusses briefly the pathway model of Mathai. 
Section 3 elaborates the closed-form representation of the thermonuclear functions for the 
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case of Maxwell-Boltzmann and a distributions in terms of G- and H- functions. Section 4 
provides conclusions. 

From Mathai and Haubold (1988) it can be seen that the expression for the reaction rate 
rij of the reacting particles i and j taking place in a nondegenerate environment is 

r l3 = n in3 {£) 2 [^y f Ea(E)e-^dE (1.1) 
= niUj(av) 

where ni and rij are the number densities of the particles i and j, the reduced mass of the 
particles is denoted by [i = J^!^. , T is the temperature, k is the Boltzmann constant, the 
reaction cross section is cr(E) and the kinetic energy of the particles in the center of mass 
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system is E = where v is the relative velocity of the interacting particles i and j. 

The reaction probability is written in the form {av ) to indicate that it is an appropriate 
average of the product of the reaction cross section and relative velocity of the interacting 
particles. For detailed physical interpretations see Haubold and Mathai (1984). 



1.1 Standard non-resonant thermonuclear function 

For non-resonant nuclear reactions between two nuclei of charges Z{ and Zj colliding at 
low energies below the Coulomb barrier, the reaction cross section has the form (Mathai and 
Haubold, 2002) 

S ( E ) -2*y(E) 

E 

where rj(E) is the Sommcrfeld parameter, h is Planck's quantum of action, e is the quantum 
of electric charge, the cross section factor S{E) is often found to be a constant or a slowly 
varying function of energy over a limited range of energy (Mathai and Haubold, 1988). The 
cross section factor S(E) may be expressed in terms of the power series expansion, 

where* = ]§and&=(^^. 



a(E) = 

with 

V (E) = 



Then 

{av) = 
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The collision probability integral, called thermonuclear function, for non-resonant ther- 
monuclear reactions in the Maxwell-Boltzmannian case is (Haubold and Mathai, 1984) 

i [-oo _ 1 

h{v,l,b,-)= / x v e~ x ~ bx "dx. (1.3) 
^ Jo 

We will consider here the general integral 

/•oo 

ii(7- l,a,b,p) = / x 7-1 e -aa; " te "' , da;, a > 0, 6 > 0, 7 > 0, p > 0. (1.4) 
Jo 



1.2 Non-resonant thermonuclear function with high energy 
cut-off 

Usually, the thermonuclear fusion plasma is assumed to be in thermodynamic equilibrium. 
But if there appears a cut-off of the high energy tail of the Maxwell-Boltzmann distribution 
function in (1.3) then the thermonuclear function is given by 

I^{y,\,b,\) = I x v e~ x - bx ~ h dx, b > 0,7 > 0,1/ > 0. (1.5) 
Again we consider the general form of the integral (1.5) as 

4 d) (7- l,a,b,p) = / x r - 1 e- ax - bx "dx, a > 0, b > 0, 7 > 0, p > 0. (1.6) 
J 

For physical reasons for the cut-off modification of the Maxwell-Boltzmann distribution func- 
tion of the relative kinetic energy of the reacting particles refer to the paper Haubold and 
Mathai (1984). 

1.3 Non-resonant thermonuclear function with depleted tail 

A depletion of the high energy tail of the Maxwell-Boltzmann distribution function of the 
relative kinetic energies of the nuclei in the fusion plasma is discussed in Haubold and Mathai 
(1984). 

For the thermonuclear function, in comparison to the strict Maxwell-Boltzmannian case, 
we have the integral 

1 f'°° x - 1 

h{y,\,l,8,b,-) = J x»e- x - x ~ bx 5 dx, 6>0, 7 >0. (1.7) 

We will consider the general integral of the type 

I 3 (j-l,a,z,6,b,p) = / sT-V *-** bx "dx, (1.8) 

Jo 
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where z > 0, a > 0, b > 0, 7 > 0, p > 0. 

1.4 Non-resonant thermonuclear function with screening 

The plasma correction to the fusion process due to a static or dynamic potential, ie the 
electron screening effects for the reacting particles, the collision probability integral to be 
evaluated in the case of the screened non-resonant nuclear reaction rate is (see Haubold and 
Mathai, 1984) 

/•OO 1 

h(is,l,b,t,p) = / x v e- x - b( - x+t ^6x. (1.9) 
Jo 

In this case we consider the general integral as 

/•OO 1 

h(j- l,a,b,t,p) = / ^"V^-^+^da;, t > 0, p > 0, 5 > 0, a > (1.10) 
Jo 

where t is the electron screening parameter. 

In the following we are evaluating the thermonuclear reaction probability integrals by a 
using pathway model (Mathai, 2005). In Section 2 we give the basic definitions that we re- 
quire in this paper. We evaluate the integrals I\ and I2 by implementing Mathai's pathway 
model and represent each of them in terms of H- function and G-function in Section 3. 

2 Mathematical Preliminaries 



We need some basic quantities for our discussion, which will be defined here. The gamma 
function denoted by T(z) for complex number z is defined as 

/•oo 

T(z)= t z ~V*dt,K(z) > (2.1) 



where K(-) denotes the real part of (•). In general T(z) exists for all values of z, positive or 
negative, except at the points z = 0, —1, —2, • • • . These are the poles of T(z). But the integral 
representation holds for the real part of z to be positive. Another important result we need 
is the multiplication formula. If z is any complex number, z 7^ 0, —1, —2, ... and let m be a 
positive integer then the multiplication formula for gamma functions is 

Vimz) = (2^) i ^m m2 -^r(z)r ( z + -)---t(z + . (2.2) 

\ m J \ m J 

For m = 2 we get the duplication formula for gamma functions, 

T(2z) = Tr-h 2z - 1 T{z)T (^z + ^j . (2.3) 
The Mellin transform of a real scalar function f(x) with parameter s is defined as 

/•OO 

M f (s) = / x s ~ 1 f(x)dx (2.4) 
Jo 
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whenever Mf{s) exists. If fi(x) and fziy) are integrable functions on the positive reals and if 
x k fi(x) and y k f 2 (y) are absolutely integrable, then the Mellin convolution property is defined 

as 

h(u) = / -h(x)f 2 (-)dx, (2.5) 
Jo X X 

then the Mellin transform of f% denoted by Mf 3 (s) is 

M f ,{s) = M h {s)M h {s) (2.6) 
where M fl (s) = f °° X s ' 1 h(x)dx and M h (s) = f °° y^ 1 f 2 (y)dy. 

The type-1 beta integral is defined as 

[ - xf-'dx = [ /-(! - yf-^y = ™, (2.7) 

R(a) > 0,K(/3) > 0. 
The type-2 beta integral is defined as 

r x a '\l + xy^dx = r y?-\l + y )-(«+«dx = ^^r, (2.8) 
Jo Jo I (a + 

»(a) > 0,»(/?) > 0. 

The G-function which is originally due to C. S. Meijer in 1936 (See Mathai, 1993; Mathai 
and Saxena, 1973) is defined as a Mellin-Barnes type integral as follows. 

where i = \/— T, L is a suitable contour and z ^ 0, m,n,p,q are integers, < m < g and 
< n < p, 

{HT=i T (Pi + s)\ (n-=ir(i - a, - ,)} 

= p J L , f } —^ (2.10) 

{nj= m+ i r(i - bj - s)} {nu+i r K' + »)} 

the empty product is interpreted as unity and the parameters ai , a 2 , ■ • • , a p and &i , 62 , • • • , 
are complex numbers such that no poles of T(bj + s), j = 1, • • • , m coincides with any pole 
of r(l — a/c — s), fe = 1, ■ ■ ■ , n; 

— bj — v / 1 — afc + A, j = 1, • • • , m; fe = 1, ■ ■ ■ , n; 1/, A = 0, 1, ■ ■ ■ . 

This means that a^ — bj 7^ I/+A+1 or cik — bj is not a positive integer for j = 1, • • • ,m; fe = 
1, ■ ■ ■ , n. We also require that there is a strip in the complex s-plane which separates the poles 
of T(bj + s), j = 1, ■ ■ ■ , m from those of T(l — — s), k = 1, • • • , n (See Mathai (1993) for 
the conditions of validity). The H- function is defined as 

Trm,n ( 1(11,0:1), ••• ,(a v ,Oi p )\ _ tt ( \ 
n P,Q \ Z \(b 1 ,f3 1 ),-,(b q ,p q ) ) ~ n \ Z > 

1 f {n^im+Mlin^ira-a.-M} 

= _1_ / L . LA . I z ~ s ds (2 11) 

2m Jc {nl= m+ i r(l - b 3 - jS )} [U% n+ i + ay*)} 



which is a Mellin-Barnes type integral where < m < q, < n < p, ctj (j = 1, 2, ■ ■ ■ ,p), f3j (j = 
1, 2, ■ ■ ■ , g) are positive real numbers and a,j (j = 1, 2, ■ ■ ■ , p) and 6j (j = 1, 2, • • • , g) are com- 
plex numbers. Any empty product is considered as unity and it is assumed that no poles of 
F(bj + (3js) for j = 1, 2, ■ ■ ■ , m coincides with any pole of T(l — aj — ctjs) for j = 1, 2, • • • , n. 
Furthermore, C is a contour in the complex s-plane moving from c—ioo to c+ioo for some real 
number c such that the points s = -4*-^ for j = 1, 2, ■ ■ ■ , m and i/ = 0, 1, ■ ■ ■ and the points 

s = 1- ^ +A for k = 1, 2, ■ ■ ■ , n and A = 0, 1, • • • lie to the left and right of C, respectively. For 
details see Springer (1979) and Mathai (1993). 

Next we utilize the pathway model of Mathai (2005). When fitting a model to experimental 
data very often one needs a model for a distribution function from a given parametric family, 
or sometimes we may have a situation of the right tail cut-off. In order to take care of 
these situations and going from one functional form to another, a pathway parameter a is 
introduced, see Mathai (2005) and Mathai and Haubold (2007). By this model we can proceed 
from a generalized type-1 beta model to a generalized type-2 beta model to a generalized 
gamma model when the variable is restricted to be positive. For the real scalar case the 
pathway model is the following, 

f(x) = c|x| 7_1 [l - o(l - a)\x\ s ]^, a > 0,5 > 0,1 - o(l - a)\x\ s > 0,7 > 0, (2.12) 

where c is the normalizing constant and a is the pathway parameter. When a < 1 the model 
becomes a generalized type-1 beta model in the real case. This is a model with the right tail 
cut-off. When a > 1 we have 1 — a = —(a — 1), a > 1 so that 

f(x) = c|x| 7_1 [l + a(a - l)\x\ 5 r^, (2.13) 

which is a generalized type-2 beta model for real x. When a — > 1 the above two forms will 
reduce to 

f(x) = c\xp~ 1 e- axS . (2.14) 
Observe that the normalizing constant c appearing in (2.12), (2.13) and (2.14), are different. 



3 Evaluation of the Integrals of Thermonuclear 

Functions 

Now we will evaluate the integrals Iiand I2 by introducing the pathway model. 

3.1 Evaluation of I\ 

r<00 

y— 1 ax— bx~ p 



e -ax-bx - dX) a > 0, 6 > 0, 7 > 0, p > 0. 



Replace e ax by [1 + a(a — l)x] As a — ► 1, [1 + a(a — l)x] becomes e ax so 

that we can extend the integral I\ to a wide class through the pathway parameter a. Let us 
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denote the wider class of this integral by I± a . 

f°° 1 

I la = ar 7-1 [l + a(a - l)z]~^e _ta 'cLc 
Jo 



where a > 1, a > 0, 5 = 1, 1 + a(a — l)x > 0, p > 0, b > 0. This is the product of two 
integrable functions. Hence we can apply Mellin convolution property for finding the value of 
the integral. Here let us take 

r / \ ) x 7 [l + aia — l)x]~ for < x < oo, a > 0, 7 > 0, a > 1 to -w 

h{x) = < 1 v 7 J - ' (3.1) 

I L), otherwise. 

, , _ f e"f P for < y < 00, p > 
/2(y) " \ 0, otherwise. (3 " 2) 

From (2.5) we have, 

r°° 1 11 

ha = / -/l(x)/ 2 (-)dx 

Jo X X 

poo 1 

= / x~ / - 1 [l + a(a-l)xr — e ~ uPx ~ P dx 
Jo 

POO 1 

= / x^' 1 [I + a(a - l)x] T^T) e - bx P dx wheie b = u p , a > I. (3.3) 
Jo 

The Mellin transform of = I\ a is the product of the Mellin transforms of f\ and / 2 

Mfy(s) = M h {s)M f2 {s) 

f°° , i_ 

M fl {s) = / z s_ V[l + a(a - l)x] («-i)dx 
JO 

/•OO ! 

= / x 7+s ^ 1 [l + a(a-l)x]"^ T )dx. 
■/ 

Putting a(a — l)x = i, we get 

1 f 00 1 

/lW [a(a- l)] 7+s 7 

1 r( 7 + *)r(^-7-*) 



[a(a - l)] 7+s 



(-) 



(3.4) 



where ^(7 + s) > 0, ^(^rr - 7 - s) > 0, a > 1. 

M h (s)= r y s ~ l ^y p dy. 
Jo 

Putting y p = t,we get 

1 f'^ 

M h {s) = - I t£ _ Vdt 
P Jo 



From (3.4) and (3.5) 

M '- w - p [«( a -i)]7+. 

Then the density of u denoted by fs(u) is available from the inverse Mellin transform. 

h(u) = 1 7 - r- ^r r +l0 °r( 7 + s )r (^__ 7 _ S W-) 

p[a(a - 1)]TT f^M 2m Jc-ioo \a-l ) \p ) 

x [a(a - l)u]^ds. (3.6) 

Comparing (3.3) and (3.6) 



r°° i 

/ x 7_1 [l+a(a-l)xpi=Te- te 'dx 
./o 

= 1 , , Hl>l (a(a - 1)^L (1 ^ +7,1) 



, , a > 1. 



Therefore, 



T 1 rr2,l f / r \?, i |( 1 "^ +7 ' 1 ) 

= 7 \H, a (a - , i un i \ 



where a > 0, b > 0, 7 > 0, p > 0, a > 1, 5R(s) > 0, 5ft(7 + s) > 0. Note that when a -> 1, /iq, 
becomes ii. But Ji Q contains all neighborhood solutions for various values of a for a > 1. 

Special Case: 

If i is an integer then put - p = m. Using equation (2.2) we get, 

^(27r)^ / Q (q-1)&™ l--l- +7 

lQ K « - D]-r (^) 1>m+1 V - m "o--.^ 

where a>0, 6 > 0, 7>0, a > 1. 

In the thermonuclear function for a non-resonant thermonuclear reaction in the Maxwell- 
Boltzmannian case 7 — 1 = 1/, a = 1, p = ^, then by using (2.3) we get, 

= [(a - !)] 7r (nil) 3 ^ 4 
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where b > 0, v > 0, a > 1 



3.2 Evaluation of I2 

1 1 

Replace e~ ax by [1 — a(l — a)x\ 1 ~ a . As a — ► 1, [1 — a(l — a)x] i-" becomes e _a:r . Let us 
denote the two integrals by and i^a respectively. 

4 d) = / xT r - 1 e - as - fce P dx, a > 0,b > 0,7 > 0,p > 0. 
■/ 

= / x 7 - 1 [l-a(l-a)x]T^e- 6 ^ P dx 

where d < a( - 1 ^ a - ) , a < 1, a > 0, 5 = 1, 1 — a(l — a)x > 0, p > 0, 6 > 0. For convenience of 

integration let us assume that d = a ^_ a ^ ■ Then is the product of two integrable functions. 
Hence we can apply Mellin convolution property for finding the value of the integral. Here 
let us take 



= Ui[l-o(l-a)i]i- for < x < ^^,0 > 0, 7 > 0, a < 1 
I 0, otherwise. 



(3.7) 



/ e~y p for < y < 00, p > 
/2(y) " \ 0, otherwise. (3 " 8) 

From (2.5) we have, 

4? = rih(x)f 2 &dx 

Jq X X 

rd 1 _ ^ 

= / x 7 " 1 [1 - a(l — a)x] < - 1 ~ a '> e~ uPx P dx where d = — - 

Jo a (l - «) 

= / ^[l-^l-aj^e^ P dx where 5 = «' , ,a<l. (3.9) 
Jo 

The Mellin transform of = is the product of the Mellin transforms of f\ and ^2- 
M f3 {s) = M h {s)M f2 {s) 

M h {s) = [ a{1 ~ a) x^ f+s ~ 1 [l-a(l-a)x]T^)dx 
Jo 

Putting a(l — a)x = i, we get 



M h {s) 



[a(l 



1 ^7 + ^(^ + 1 

Ki-a)]^ r ( 1 + 7 + _i_ + s ) 



, R( 7 + s) >0,a < 1. (3.10) 
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From equation (3.5) we get 

M /2 ( S ) = lrQ, $t(s)>0. (3.11) 

From (3.10) and (3.11) 

M r.i 1 rh + s)r (^ + 1 ) r 00 
>' {s > „Ki-«)]-r+. r ( 1+7 + _L_ + s ) ■ 

Then the density of it denoted by fs(u) is available from the inverse Mellin transform. 

/ 3 (u) = -f- — / — f4o(l - a)n]-d S . (3.12) 

p[a(l - a)]T 2vrt 7 c _ ioo r / 1+7+ _l_ + s j 

Comparing (3.9) and (3.12) 

[ d X ^- l [l - o(l - a)x]i^e- 6 ^ P dx 
JO 

1=5 + 1) 1 r( 7 + s )r^ 



/ — \t!—^[a{\ -a)b-p]- s ds 

Jc-ioo r 11 -I- 'v -I- -J- 1- .si 



p[a(l - a)p 2m J c _ ioo r ^ + 7 + _1_ + ^ 
p[a(l - a)]" 



xf(«(i-^i ( ( ;i; < ;ir' 1) )- a < 1 ' 



Therefore, 



^ " p[o(l - a)]T ^ P j J 



where a > 0, b > 0, 7 > 0, p > 0, a < 1, a — > 1. 
Special Case: 

If i is an integer then put - p =m. Then following through the same procedure as before 
one has 

M) = v^(2vr)^r(— + lj q / o(1 _ a)6 m x 

where a > 0, 6 > 0, 7 > 0, a < 1. 

For the thermonuclear function fornon-resonant thermonuclear reactions with high energy 
cut off 7 — 1 = v, a = 1, P=t? then we get, 



r(d) _ r ( 1 -" +1 ) ^3,0 f 



3,0 ( (1 - a)b 2 ,v+j^+2 
4 
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where b > 0, v > 0, a < 1. The importance of the above result is that I 2 a gives an extension 
of the integral to a wider class of integrals through the pathway parameter a, and their 
solutions. 



4 Conclusion 



In the field of stellar, cosmological, and controlled fusion, for example, the core of the 
Sun is considered as the gravitationally stabilized solar fusion reactor. The probability for a 
thermonuclear reaction to occur in the solar fusion plasma depends mainly on two factors. 
One of them is the velocity distribution of the particles in the plasma and is usually given 
by the Maxwell-Boltzmann distribution of Boltzmann-Gibbs statistical mechanics. The other 
factor is the particle reaction cross-section that contains the dominating quantum mechanical 
tunneling probability through a Coulomb barrier, called Gamow factor. Particle reactions in 
the hot solar fusion plasma will occur near energies where the product of velocity distribu- 
tion and tunneling probability is a maximum. The product of velocity distribution function 
and penetration factor is producing the Gamow peak. Mathematically, the Gamow peak is a 
thermonuclear function. In case of taking into consideration electron screening of reactions in 
the hot fusion plasma, the Coulomb potential may change to a Yukawa-like potential. Tak- 
ing into account correlations and long-range forces in the plasma, the Maxwell-Boltzmann 
distribution may show deviations covered by the distribution predicted by Tsallis statistics 
in terms of cut-off or depletion of the high-velocity tail of the distribution function. In this 
paper, closed-form representations have been derived for thermonuclear functions, thus for 
the Gamow peak, for Boltzmann-Gibbs and Tsallis statistics. For this purpose, a generalized 
entropy of order a and the respective distribution function, have been considered. The case 
a = 1 recovers the Maxwell-Boltzmannian case. This general case is characterized by mov- 
ing cut-off, respectively the upper integration limit of the thermonuclear function to infinity. 
The closed-form representations of thermonuclear functions are achieved by using generalized 
hypergeometric functions of Fox and Meijer or H-functions and G-functions, respectively. 
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